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ABSTRACT 


Experiments  on  pressure  distribution  and  heat  transfer  on 

a  grid  simulating  a  parachute  cloth  were  performed.  The  porosity 

of  the  grid  was  25%.  A  pressurized  wind-tunnel  served  as  the  flow 

facility.  The  average  approach  of  the  flow  was  134  feet  per  second. 

The  pressure  ratio  P  /P  was  considered  as  the  main  variable 

UP  down 

and  varied  from  around  11  to  1.  45.  A  transient  energy  balance  was 
used  to  calculate  local  heat  transfer. 


ill 


CONTENTS 


SECTION  PAGE 

I.  INTRODUCTION  1 

II.  EXPERIMENTAL  EQUIPMENT  4 

A.  Facility  4 

B.  Heating  Chamber  4 

C.  Instrumentation  5 

1.  Pressure  5 

2.  Temperature  5 

D.  The  Experimental  Model  6 

III.  EXPERIMENTAL  STUDIES  9 

A.  Uniformity  of  Approach  Flow  9 

B.  Surface  Pressure  Distribution  9 

C.  Recovery-Temperature  Runs  10 

D.  Heat  Transfer  Runs  13 

IV.  CONCLUSIONS  16 

REFERENCES  17 

APPENDIX  -  Transient  Energy  Balance  49 

iv 


ILLUSTRATIONS 


FIGURE  PAGE 

1  Heating  apparatus,  pneumatic  piston,  steel 

windows  for  Schlieren  pictures.  18 

2  Heating  chamber,  with  the  test  model.  19 

3  Typical  photograph  of  the  manometer  board 

used  to  read  the  pressure  distribution.  20 

4  View  of  the  test  model;  (a)  without  supporting 

pieces,  and  (b)  with  supporting  pieces.  21 

5  Detailed  sketch  of  the  test  model.  22 

6  Cross-sections  where  measurements  were 

taken.  23 

7  Typical  instrumentation  of  two  cross-sections 

corresponding  to  Bars  2  and  1,  respectively.  23 

8  Sketch  of  the  test  section  of  the  wind  tunnel.  24 

9  Pressure  distribution  for  Cross-Section  A-A.  25 

10  Pressure  distribution  for  Cross-Section  B-B.  26 

11  Pressure  distribution  for  Cross-Section  C-C.  27 

12  Pressure  distribution  for  Cross-Section  A'-A*.  28 

13  Pressure  distribution  for  Cross-Section  B'-B'.  29 

14  Pressure  distribution  for  Cross-Section  C'-C'.  30 

15  Recovery  temperature  distribution  for  Cross- 

Section  A-A.  31 

16  Recovery  temperature  distribution  for  Cross- 

Section  B-B.  32 


V 


✓ 

1 


ILLUSTRATIONS  (cont'd) 


FIGURE 

PAGE 

17 

Recovery  temperature  distribution  for  Cross- 
Section  C-C. 

33 

18 

Recovery  temperature  distribution  for  Cross- 
Section  A'-A\ 

34 

19 

Recovery  temperature  distribution  for  Cross- 
Section  B'-B'. 

35 

20 

Recovery  temperature  distribution  for  Cross- 
Section  C'-C'. 

36 

21 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  A-A. 

37 

22 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  B-B. 

38 

23 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  C-C. 

39 

24 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  A'-A'. 

40 

25 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  B'-B'. 

41 

26 

Distribution  of  the  heat  transfer  coefficients  for 
Cross-Section  C'-C'. 

42 

27 

Distribution  of  the  heat  transfer  coefficients  in 
polar  form  for  Cross-Section  A-A. 

43 

28 

Distribution  of  the  heat  transfer  coefficients  in 
polar  form  for  Cross-Section  B-B. 

44 

29 

Distribution  of  the  heat  transfer  coefficients  in 
polar  form  for  Cross-Section  C-C. 

45 

30 

Distribution  of  the  heat  transfer  coefficients  in 
polar  form  for  Cross-Section  A'-A'. 

46 

31 

Distribution  of  the  heat  transfer  coefficients  in 
polar  form  for  Cross-Section  B'-B'. 

47 

32 

Distribution  of  the  heat  transfer  coefficients  m 
polar  form  for  Cross-Section  C'-C'. 

Vi 

48 

I 

t 


ILLUSTRATIONS  (cont'd) 


FIGURE 

33 

34 


PAGE 

Typical  extrapolation  technique  for  heat 

transfer  coefficient.  56 

Typical  temperature-time  history.  57 


vii 


H  \ 


✓ 

< 


A 

a 

c. 

D 

h 

k 

P 

q 

R 

r 

T 

t 

e 

e 

or 

r 


NOMENCLATURE* 

2 

convective  heat  transfer  area,  ft. 
speed  of  sound,  ft.  / sec. 
specific  heat,  Btu/lb.  °F. 
outside  diameter  of  the  bar,  ft. 

convective  heat  transfer  coefficient,  Btu/hr.ft.  °F. 
thermal  conductivity,  Btu/hr.  ft.  °F. 
absolute  static  pressure,  inches  of  Hq. 

2 

rate  of  heat  transfer  per  unit  area,  Btu/hr.  -ft. 

dimensionless  temperature  ratio 

ratio  adiabatic  wall  temperature  to  stagnation 
temperature 

temperature,  °R 

thickness,  ft. 

radiations'  interchange  factor 

angle  measured  around  the  model  from  upstream 
direction 

viscosity,  lb. /sec.  ft. 
density,  lbm/ft3 

„  2  2 
Stephan-Boltzmann  constant,  Btu/ft.  hr.  °R 

time,  hr. 


♦  Several  symbols  defined  in  the  text  and  used  but  briefly,  are  not 
included  in  this  list. 
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NOMENCLATURE  (cont'd) 


Superscripts 

*  denotes  same  condition  calculated  from  the 

stagnation  conditions  upstream,  assuming  an 
isentropic  expansion  to  sonic  velocity 

Subscripts 

aver.  average  value 

aw  denotes  adiabatic  wall  condition 

down  denotes  conditions  downstream 

m  refers  to  model  material 

o  denotes  stagnation  conditions  in  main  stream 

stag  denotes  conditior  s  at  the  stagnation  point 

up  denotes  conditions  upstream 


w 


denotes  local  conditions  at  the  model  surface 


I.  INTRODUCTION 


This  report  describes  the  methods  and  results  of  an  experi¬ 
mental  program  to  determine  the  convective  heat  transfer  to  the  cloth 
of  a  parachute  designed  to  decelerate  vehicles  from  supersonic  speeds. 
Since  instrumentation  difficulties  precluded  the  possibility  of  deter¬ 
mining  the  distribution  of  heat  flux  to  the  fiber  of  which  the  cloth  is 
composed,  the  experimental  program  was  divided  into  two  parts;  the 
microscale  study  and  the  macroscale  study.  In  the  microscale  study, 
reported  in  this  document,  an  enlarged  model  of  the  woven  fiber  was 
built  and  instrumented  to  determine  how  the  heat  flux  is  distributed 
to  the  fiber.  In  this  model,  the  simulated  fibers  were  one  inch  in 
diameter.  In  the  macroscale  study,  to  be  reported  later,  the  gross 
convective  heat  transfer  per  unit  area  to  an  actual  scale  model  of 
the  woven  cloth  was  determined.  The  porosity  of  the  microscale-study 
grid  was  fixed  by  mass  flow  considerations  for  existing  tunnel 
facilities.  The  incoming  feed  lines  along  with  the  quick  opening  valve 
fixed  the  incoming  mass  flow.  Then  the  model  porosity  was  computed 
to  produce  a  sonic  throat  in  and  at  the  open  area  minimums.  On  this 
basis,  the  grid  was  built  with  a  porosity  of  25%.  The  macroscale- study 
grid  duplicates  the  actual  size  and  porosity  of  a  representative  parachute 


mesh  cloth. 
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The  experimental  facility  used  in  the  microscale  study  is 
located  at  the  University  of  Minnesota's  Rosemount  Research  Center. 

It  was  used  before  to  obtain  corresponding  results  reported  in  Reference 
1  on  a  model  of  a  ribbon  parachute.  For  the  present  microscale  study, 
two  differences  in  procedure  from  those  used  in  Reference  lwere  made. 
The  first  involves  the  operating  pressure  level.  In  Reference  1  the  air 
entered  the  wind  tunnel  from  a  compressor,  where  it  was  dried  to  a  dew 
point  of  -40°F,  and  was  exhausted  into  Uie  atmosphere.  In  the  present 
study,  however,  due  to  the  fact  that  the  model  is  larger  than  full  scale 
the  sonic  Reynolds  number  is  increased  proportionally.  The  sonic 


Reynolds  number  may  be  defined  si 


*  * 
r  cl  o 

- ~F~ 


Since  the  heat  transfer  depends  on  the  Reynolds  number,  essentially 
constant  throughout  this  study,  it  was  desirable  to  compensate  for 
this  increase  by  reducing  the  pressure  and  hence  the  density  of  the  flow. 

For  this  reason,  it  was  decided  to  take  air  directly  from  the  atmosphere 
and  exhaust  it  into  a  vacuum  tank.  With  this  change,  the  Reynolds  number 
was  between  700,  000  and  900,  000.  Unfortunately,  in  doing  this  it  was 
inconvenient  to  dry  the  incoming  air.  In  the  first  experimental  runs  made  to 
determine  the  heating  rate  of  the  model,  the  necessary  difference  between  the 
initial  model  temperature  and  the  steady  state  (recovery)  temperature  was 
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obtained  in  the  same  way  as  in  Reference  1.  In  this  method  the  model 
was  cooled  with  dry  snow  (CO^)  to  a  temperature  of  about  -110°F  before 
inserting  it  into  the  wind  tunnel;  however,  in  the  present  experiments  the 
moisture  in  the  undried  air  was  often  seen  to  produce  a  frost  layer  on  the 
model.  Since  this  frost  layer  would  influence  the  results  by  an  unknown 
amount,  it  was  decided  to  revise  the  procedure  (the  second  difference  from 
Reference  1)  by  heating  the  model  prior  to  insertion  into  the  wind  tunnel. 
This  procedure  was  actually  found  to  be  quicker  and  still  provided  a 
sufficiently  uniform  initial  temperature. 

In  the  following  sections,  the  important  pieces  of  special  equip¬ 
ment  used  in  the  tests  and  the  model  itself  are  described  first.  Following 
this,  measurements  of  the  pressure  distribution  around  the  model  are 
discussed.  Then  the  two  main  parts  of  the  experimental  program  are 
described  separately.  First  are  the  recovery-temperature  measurements 
needed  to  find  the  steady-state  temperature  of  the  model  in  the  flowing 
stream,  and  second  are  the  transient  temperature  measurements  needed 
to  find  the  heating  rate.  The  details  of  the  procedure  used  to  compute  the 
heat  transfer  coefficients  from  the  data  of  the  two  types  of  measurements 
are  given  in  the  Appendix.  Data  were  taken  at  six  different  pressure 
ratios.  These  varied  from  eleven,  which  was  the  highest  obtainable 
with  the  facility  used,  to  1.45,  which  was  as  close  to  unity  as  was  felt  to 


be  of  interest. 
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II.  EXPERIMENTAL  EQUIPMENT 

A.  Facility 

The  experimental  facility  used  is  the  same  as  described  in  Section 
ILA  of  Reference  1,  and  shown  in  Figure  2  of  that  reference  with  the 
following  exceptions: 

1.  Upstream  of  the  teat  section  the  air  was  taken  directly 
from  the  atmosphere  and  therefore  the  compressed  air 
system  was  not  used. 

2.  Downstream  of  the  model  air  was  led  into  the  vacuum  tanks. 
The  vacuum  system  has  a  capacity  of  22,750  cubic  feet 

and  is  evacuated  by  two  Allis-Chalmers  27-D  vacuum 
pumps.  These  pumps,  rated  at  3390  cfm  each,  are  powered 
by  a  350  h.  p.  motor,  which  drives  both  pumps.  The  lowest 
vacuum  obtainable  with  the  pumps  in  series  was  2"  Hg. 

B.  Heating  Chamber 

To  carry  out  the  heat  transfer  measurements  it  was  necessary  to 
heat  the  model  to  a  constant  temperature  before  inserting  it  into  the  air 
stream.  Figures  1  and  2  show  two  different  views  of  the  apparatus 
constructed  to  do  this.  The  pneumatic  piston  in  these  figures  made 
possible  the  insertion  and  removal  of  the  model  as  fast  as  necessary.  By 
regulating  the  valves  it  was  possible  to  avoid  destructive  impacts  against 


the  wind-tunnel  windows.  The  chamber  exterior  was  built  with  removable 
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top  and  bottom  plates  to  facilitate  handling  of  the  heater.  The  heater 
consists  of  coiled  Nichrome  wire  uniformly  distributed  on  two  removable 
asbestos  plates  placed  on  the  side  walls  of  the  heating  chamber.  During 
the  heating  period  the  top  plate  was  kept  closed.  The  Dymec  system 
described  in  Subsection  C,  monitored  the  temperature  distribution  during 
this  period. 

C.  Instrumentation 

1.  Pressure 

In  order  to  measure  the  pressure  distribution  the  mercury 
manometer  board  shown  in  Figure  2  of  Reference  1  was  used. 
Pressure  differences  could  be  read  with  an  uncertainty  of  +  0.025 
inches  of  mercury.  A  quick  acting  valve  at  the  top  clamped  the 
leads  from  the  model,  thus  holding  the  reduced  pressure.  At  this 
moment  readings  were  taken  by  photographing  the  manometer 
board  so  that  the  readings  could  be  rechecked.  One  of  these  photo¬ 
graphs  is  shown  in  Figure  3.  The  pressure  taps  were  silver- 
soldered  in  holes  0.  028  inches  in  diameter  and  drilled  in  the  wall. 
Figure  7  gives  a  detailed  explanation.  The  outside  surface  was 
repolished  after  this  installation.  All  pressure  taps  were  tested 
and  found  to  be  leakproof. 

2.  Temperature 

Thirty  gauge,  roll  calibrated  iron-constantant  thermocouples 
were  inserted  approximately  i/32  inch  in  the  wall  of  the  model. 
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After  this  the  thermocouples  were  silver-soldered. 

Modification  M24  of  the  Dymec  DY-2010A  Data  acquisitor's 
system  provided  25  channel  inputs  continuously  monitored  at  a 
precise  rate  of  5  channels  a  second,  with  digital  read  out. 

D.  The  Experimental  Model 

Two  general  views  of  the  test  model  or  grid  are  shown  in  Figure  4. 
It  consists  of  seven  highly  polished  bars  of  free  machining  stainless 
steel,  arranged  in  such  a  way  that  they  resemble  woven  cloth  of  25% 
porosity.  Stainless  steel  was  chosen  because  of  its  relatively'  small 
thermal  conductivity  and  corrosion  resistance.  It  should  be  recalled, 
however,  that  since  only  the  rate  of  heat  transfer  through  the  air  to  the 
surface  is  being  considered,  the  material  used  in  the  model  can  be  chosen 
at  the  discretion  of  the  experimenter. 

Bars  1,  2,  and  3,  Figure  5,  are  hollow  with  1/16"  wall  thickness. 
Bar  1  contains  the  pressure  taps  and  Bar  2  the  thermocouples.  Because 
of  the  influence  of  the  thickness  of  these  hollow  bars  in  the  calculation 
of  the  heat  transfer  coefficients,  the  bars  were  made  as  accurate  as 
possible  within  a  tolerance  of  +  .  003  inches.  Bar  3  was  also  made 
hollow  because  it  is  in  contact  with  Bars  1  and  2,  and  the  rest  of  the  bars 
are  solid.  Bars  1,  2,  and  3  were  assembled  by  cementing  together  two 
parts  which  were  mirror  images. 


The  model  is  held  together  by  two  thin  plates,  one  at  the  top 
and  one  at  the  bottom  (see  Figure  4b),  and  two  heavier  pieces  which  fit 
snugly  into  steel  "windows"  in  the  two  opposite  wind-tunnel  walls  as 
sho'v n  in  Figure  1.  With  this  construction,  it  was  possible  to  retract  the 
model  from  the  wind-tunnel  to  heat  it  up  to  the  desired  temperature.  The 
optical  windows  shown  on  Figure  1  are  provided  for  shadowgraph/ schlieren 
observations  of  the  flow. 

With  respect  to  the  orientation  of  the  model  with  the  direction  of 
flow,  measurements  were  taken  at  the  six  circular  cross  sections  indicated 
in  Figure  6.  The  necessary  information  was  obtained  from  five  of  these 
cross  sections  while  the  sixth,  C-C  or  C'-C,  served  to  check  the  symmetry 
of  the  flow.  This  symmetry  is  more  apparent  from  Figure  5  which  shows  the 
proximity  of  the  instrumented  bar  with  respect  to  all  neighboring  perpen¬ 
dicular  bars.  As  mentioned  above,  however,  Bar  1  was  instrumented  only 
with  pressure  taps  and  Bar  2  only  with  thermocouples.  Readings  at  all 
six  cross  sections  shown  in  Figure  6  were  therefore  obtained  by  rotating 
the  model  180°  and  repeating  the  run.  Thus  cross  sections  A-A,  B-B,  and 
C-C  are  transformed  into  cross  sections  A'-A',  B'-B',  and  C'-C', 
respectively  by  turning  the  model  around. 

The  symmetry  of  the  flow  with  respect  to  the  horizontal  axes  X-X 
and  Y-Y,  respectively,  shown  in  Figure  7,  made  it  necessary  to  place 
instrumentation  around  only  half  of  the  cross  sections.  The  half  cross 
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sections  instrumented  were  those  farthest  from  the  wind-tunnel  walls. 

The  cross  sections  are  instrumented  at  seven  peripheral  locations 
corresponding  to  the  angles  0  =0°,  30°,  60°,  90°,  120°,  150°,  and  180° 
measured  from  the  upstream  side  with  the  exception  that  in  cross  sections 
A-A  and  A'-A'  there  is  no  instrumentation  at  the  point  of  contact.  See 
Figure  7. 
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III.  EXPERIMENTAL  STUDIES 
A.  Uniformity  of  Approach  Flow 

According  to  the  measurements  taken  in  the  summer  of  1964 
by  Eckert  and  Scott  on  the  aame  wind  tunnel,  the  upatream  velocity 
distribution  in  as  given  in  Figure  8  of  Reference  1.  All  velocities  are 
within  1%  of  th ;  centerline  velocity  with  the  exception  of  the  four  points 
taken  near  the  nozzle  walls. 


B.  Surface  Pressure  Distribution 

The  pressure  distribution  around  the  model  in  the  cross-sections 

A -A,  B-B,  and  C-C  was  measured  at  different  pressure  ratios 

P  /p  .  .  The  model  was  then  turned  180°  around  a  horizontal  axis 

up'  down 

through  its  mid-point,  as  explained  in  Section  IID,  and  the  runs  were 
repeated.  In  this  case,  the  pressure  distribution  was  measured  at  the 
cross  sections  A'-A',  B’-B',  and  C'-C'. 

Since  the  air  was  taken  directly  from  the  atmosphere  through  a 
piping  system,  the  total  pressure  P  varied  only  with  the  atmospheric 
conditions.  One  wall  static  pressure  tap  located  upstream,  (see  Figure  8), 


where  the  velocity  of  the  air  is  negligible,  gave  P 
wall  static  pressure  tap  located  in  the  centerline  of  the  last  window  (see 


up 


p  .  Another 
stag 


Figure  8),  was  used  to  measure  Pj0wil» 
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Figures  9  though  14  ahow  the  ratio  P/P  aia  function  of  the 

up 

angle  9,  for  the  aix  croaa-aectiona.  The  almost  conatant  pressure 

distribution  around  the  cross- section  A -A,  Figure  9,  can  be  explained  by 

the  fact  that  the  flow  is  retarded  by  the  presence  of  the  rearward  bar, 

thus  causing  the  pressure  to  remain  near  the  stagnation  pressure.  The 

rest  of  the  figures  show  that  the  pressure  distribution  around  the  upstream 

portion  of  the  cross-section  is  independent  of  pressure  ratio  up  to  a 

certain  value  of  the  angle  9,  say  9^.  This  can  be  explained  on  the  basis 

of  Schlichting's  theory  of  separation  (Reference  2).  This  theory  states 

that  the  angle  9  at  which  separation  occurs  increases  with  pressure  ratio, 

P  / P.  .  Thus  we  can  infer  that  the  point  9.,  at  which  the  pressure 
up  down  • 

distribution  curve  for  a  given  pressure  ratio  breaks  away  from  the  curves 
for  higher  pressure  ratios,  is  the  separation  point. 

Some  differences  in  the  pressure  distributions  given  by  Figures 
11  and  14,  which  correspond  tocross-sections  having  the  same  orientation 
with  respect  to  the  flow  are  probably  due  to  the  non-perfect  symmetry  of 
the  flow  and  to  the  non-perfect  alignment  of  the  model  bars. 

C.  Recovery  Temperature  Runs 

To  calculate  the  heat  transfer  coefficients  by  the  method 
described  in  the  Appendix,  the  recovery  temperature,  also  called 
adiabatic  wall  temperature,  must  be  known.  The  recovery  temperature 
at  each  point  around  a  surface  is  the  temperature  that  this  point  assumes 


-11- 

in  steady  state  under  the  influence  of  internal  friction  in  the  boundary 
layer  where  the  convective  heat  flux  has  the  value  zero. 

At  the  highest  pressure  ratio,  pup^pdown^  the 

characteristic  running  time  of  the  wind  tunnel  was  about  4  minutes. 
After  this  time  itwaj  not  possible  to  maintain  a  constant  pressure  ratio 
since  the  pumping  system  was  not  able  to  pump  out  the  air  flowing  into 
the  vacuum  tanks.  Because  of  this  the  downstream  pressure  started  to 
increase  and  the  pressure  ratio  decreased.  The  rate  of  change  of  the 
pressure  ratio  decreased  with  decreasing  pressure  ratio. 

Data  on  recovery  temperatures  were  taken  at  five  different 
pressure  ratios.  These  data  give  sufficient  information  because  of  the 
small  changes  of  the  recovery  temperature  with  the  preosure  ratio  as 
explained  below.  Two  of  these  pressure  ratios  are  supercritical,  two 
are  subcritical  and  the  fifth  one  close  to  the  critical  pressure  ratio. 

In  taking  the  data  at  the  highest  pressure  ratio,  it  could  be 
observed  that,  after  the  4-minute  interval,  during  which  the  pressure 
ratio  remained  constant,  the  rate  of  change  of  the  temperature  was  at 
all  points  less  than  .02°F/sec.  Hence  the  last  temperatures  measured 
were  taken  as  the  recovery  temperatures  and  no  extrapolations  were 
necessary.  To  obtain  the  recovery  temperature  at  the  lower  pressure 
ratios,  the  wind  tunnel  was  kept  running  beyond  the  4-minute  interval. 
As  the  pressure  ratio  gradually  decreased,  each  set  of  temperature 
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readings  were  taken  at  different  pressure  levels  symmetrical  around 
each  desired  pressure  ratio;  the  first  one  before,  the  second  one  at, 
and  the  third  one  after  reaching  the  desired  pressure  ratio.  As  all  three 
of  these  sets  of  readings  showed  nearly  the  same  temperatures,  these 
temperatures  were  taken  as  the  recovery  temperatures.  This  meant 
essentially  that  the  time  constant  for  temperature  relaxation  was  small 
compared  to  the  time  rate  at  which  the  pressure  changed.  After  20 
minutes  running  time  the  lowest  pressure  ratio  was  reached.  By  plotting 
Taw/TQ  vs.  Pup/Pjown  **  wa*  *oun<*  t*'at  the  Poin**  corresponding  to  the 
highest  and  lowest  pressure  ratio  were  different  by  less  than  2%.  This 
supported  the  idea  of  accepting  the  temperature  measured  in  the  manner 
indicated  as  the  recovery  temperature.  A  very  slight  increase  of  the 
recovery  temperature  with  decreasing  pressure  ratio  could  be  observed. 

During  the  recovery  temperature  runs,  data  were  taken  at  three 
different  cross  sections  at  the  same  time.  A  total  of  21  thermocouples  were 
used,  20  of  them  at  the  three  cross  sections  and  the  other  one  placed 
upstream  to  measure  the  stagnation  temperature  of  the  air  T  .  This 
last  one  had  two  connections  to  the  Dymec,  so  that  it  was  possible  to 
get  more  eccurate  information  about  the  variations  of  the  air  stagnation 
temperature  during  the  4.4  seconds  required  for  the  Dymec  to  take  the 
readings. 

The  temperature  variations  of  the  air  stagnation  temperature 
have  been  less  than  +  0.  5°F,  during  each  run. 
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Each  run  was  repeated  and  showed  consistent  results.  The 
ratio  T  ^/Tq  of  the  absolute  temperatures  has  been  plotted  as  functions 
of  the  angle  0  in  Figures  15  through  20.  Cross  sections  A-A,  B-B, 
and  C-C  indicated  higher  ratios  of  T  /T  at  corresponding  angles  0, 
than  cross  sections  A'-A',  B'-B',  and  C'-C'.  Cross  sections  C-C  and 
C'-C1  should  show  the  same  results.  The  difference  may  be  attributable 
to  a  probable  change  of  the  relative  humidity  of  the  air,  to  the  non-perfect 
symmetry  of  the  flow,  and  to  non-perfect  alignment.  The  figures  indicate 
that  the  temperature  recovery  is  almost  100%  over  the  upstream  part  of 
the  model  as  well  as  in  the  separated  flow  regime. 

D.  Heat  Transfer  Runs 

Iu  order  to  get  the  desired  information  on  the  local  heat  transfer, 
transient  data  were  also  taken  in  the  cross  sections  mentioned  above.  The 
general  heat  transfer  experimental  technique  is  described  in  detail  in 
the  Appendix.  As  mentioned  above,  the  technique  required  preheating 
the  model  with  an  electric  heater  to  a  temperature  significantly  diffe  ent 
from  the  recovery  temperature.  To  avoid  circumferential  heat  conduction 
in  the  model  it  was  more  important  to  obtain  a  uniform  temperature 
distribution  at  the  start  of  each  run  than  to  get  a  prescribed  temperature. 
By  observing  the  temperature  readings  on  the  Dymec  during  the  heating 
period,  it  was  possible  to  initiate  the  run  when  the  temperatures  were 
sufficiently  uniform.  In  general  the  model  was  heated  to  around  250°F. 
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The  largest  temperature  difference  in  the  cross  section  was  always 
less  than  10°F.  The  conduction  error  in  the  heat  transfer  coefficients 
produced  by  this  difference  was  estimated  and  found  negligible.  Once 
the  temperature  distribution  was  uniform  to  this  degree,  the  model  was 
inserted  in  the  wind  tunnel  in  less  than  1/2  second.  At  the  same  time 
two  operators  opened  the  upstream  quick-opening  DC  Zurik  valve  and 
the  Butterfly  valve  downstream  between  the  wind  tunnel  and  the  vacuum 
tanks.  In  total,  the  system  was  operating  in  less  than  1  second.  Owing 
to  the  rapid  cooling  of  the  grid  upon  insertion  into  the  tunnel,  rapid 
measurements  of  wall  temperatures  were  necessary.  Thus  only  one 
cross-section  was  studied  in  each  heat  transfer  run. 

The  heat  transfer  coefficients  were  extrapolated  to  zero  time 
ignoring  data  taken  in  the  first  second,  since  this  data  gave  unreliable 
results  arising  from  disturbances  created  by  the  actual  insertion  of  the 
test  grid.  The  reader  is  again  referred  to  the  Appendix  where  a  sample 
run  and  the  various  steps  in  the  data  reduction  are  illustrated. 

Data  were  taken  at  six  different  pressure  ratios,  three  of  them 
corresponding  to  supercritical  pressure  ratios,  two  to  subcritical 
ratios,  and  the  sixth  one  near  the  critical  pressure  ratio. 

Local  heat  transfer  coefficients  are  presented  in  Figures  21 
through  26  in  Cartesian  coordinates  and  in  Figures  27  to  32  in  polar 
diagrams.  The  influence  of  the  pressure  ratio  on  heat  transfer  is 
somewhat  puzzling.  This  may  be  due  to  the  complicated  geometry  and 
also  possibly  to  the  fact  that  the  humidity  in  the  air  was  rather  high  on 


-15- 


days  at  which  some  of  the  experiments  were  run.  On  the  other  hand, 
the  pressure  distribution  curves  in  Figures  11  through  16  exhibit  quite 
a  regular  decrease  of  pressures  with  increasing  pressure  ratio  and 
show  no  indication  of  a  humidity  influence.  To  draw  curves  for  each 
pressure  ratio  in  the  figures  would  have  confused  the  picture.  Therefore 
only  an  average  curve  has  been  indicated  in  Figures  21  through  23  and  25. 
The  measured  points  in  Figures  24  and  26  separate  clearly  into  two 
groups:  one  for  subcritical  and  the  other  for  supercritical  pressure 
ratio.  Two  curves  have  therefore  been  drawn  in  these  figures. 

Figures  27  through  32  repeat  the  curves  contained  in  the  preceding 
figures  in  polar  diagrams.  One  or  two  scales  have  been  added  which 
allow  the  ratio  of  local  to  average  heat  transfer  coefficients  to  be  read. 

It  is  felt  that  this  ratio  is  especially  useful  in  design  calculations.  The 
value  of  the  average  heat  transfer  coefficients  itself  can  easily  be  read 


from  the  left-hand  scale. 
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IV.  CONCLUSIONS 

Experiment!  have  been  performed  to  determine  local 

convective  heat  transfer  to  a  model  simulating  the  shape  of  parachute 

cloth  with  a  porosity  of  25%.  The  pressure  ratio  P  /P  was 

UP  down 

varied  from  subcritical  to  supercritical  values  in  the  range  of  1.45  to 
approximately  11.  The  upstream  total  pressure  was  close  to 
atmospheric. 

Pressure  distribution,  recovery  temperatures,  and  convective 
heat  transfer  coefficients  have  been  obtained  and  are  presented  in 
Figures  9  through  32.  The  complicated  nature  of  the  flow  prevented 
any  analytical  treatment. 
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FIGURE  1  KEATING  APPARATUS.  PNEUMATIC  PISTON.  STEEL  WINDOWS  TOR 
^HLIEREN  PICTURES 


FIGURE  3  TYPICAL  PHOTOGRAPH  OF  THE  MANOMETER  BOARD  USED  TO  READ  TOE  PRESSURE 


FIGURE  4a  VIEW  OF  THE  TEST  MODEL  WITHOUT 
SUPPORTING  PIECES 


FIGURE  4b  VIEW  OF  THE  TEST  MODEL  WITH  SUPPORTING 
PIECES 
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FIGURE  5  DETAILED  SKETCH  OF  THE  TES"1  MODEL 


FIGURE  8  SKETCH  OF  IHE  TEST  SECTION  OF  TEE  WIND  TUNNEL 


FIGURE  9  PRESSURE  DISTRIBUTION  FOR  CRQflB  SECTION:  A -A 


FIGURE  13  PRESSURE  DISTRIBUTION  FCR  CROSS-SECTION  B*-B* 


FIGURE  1$  RECOVER!  TEMPERATURE  DISTRIBUTION  FCR  CROSS  ^SECTION  A -A 
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FIGURE  16  REC07EHX  TEMPERATURE  DISTRIBUTION  FOR  CROSS-SECTION  B-B 


001 


FIGURE  18  HEC07EHT  TEMPERA TUIE  DISTRIBUTION  FOR  CROSS-SECTION  A»-A» 


FIGURE  20  RECOVER!  TEMPERATURE  DISTRIBUTION  t’GR  CROSS -SECTION  O' -O' 
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FIGURE  21  DISTRIBUTION  CF  THE  HEAT  TRANSFER  COEFFICIENTS  FOR  CROSS-SECTION  JU 
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FIGURE  23  DISTRIBUTION  OF  THE  HEAT  TRANSFER  COEFFICIENTS  FOR  CROSS -SECTION  C-C 
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FIGURE  2 U  DISTRIBUTION  GF  THE  HEAT  TRANSFER  COEFFICIENTS  FOR  CROSS-SECTION  A*  Ji* 


FIGURE  25  DISTRIBUTION  OF  THE  HEAT  TRANSFER  COEFFICIENTS  FOR  CROSS-SECTION  B'-B' 
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Cro$8-Sect?on  A- A 
90*  • 


FIGURE  27  DISTRIBUTION  OF  THE  HEAT  TRANSFER  COEFFICIENTS  IN  TOLAR  FORM 
FOR  CROSS-SECTION  A -A 
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Crost- Section  B-B 


FIGURE  28  DISTRIBUTION  OF  THE  HEAT  TRANSFER  COEFFICIENTS 
IN  POLAR  FORM  FOR  CROSS-SECTION  B-B 
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FIGURE  29  DISTRIBUTION  CF  THE  HEAT  TRANSFER  COEFFICIENTS 
IN  POLAR  FORM  FOR  CROSS-SECTION  CX 
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Figure  30  distribution  of  the  heat  transfer  coefficients 

IN  POLAR  FORM  FOR  CROSS-SECTION  A'^A' 
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Cross -Section  B'-B* 


FIGURE  31  DISTRIBUTION  CF  THE  HEAT  TRANSFER  COEFFICIENTS  IN 
POLAR  FORM  FCu  CROSS -SECTION  B'-B* 
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FIGURE  32  DISTRIBUTION  OF  THE  HEAT  TRANSFER  COEFFICIENTS 
IN  POLAR  FORT  FOR  CROSS-SECTION  C'-C' 
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APPENDIX 

TRANSIENT  ENERGY  BALANCE 
1.  Basic  Equation 

The  transient  technique  relies  on  the  heat  capacity  of  the  model 
skin  for  the  experimental  determination  of  convective  heat  flux.  The 
general  heat  balance  on  a  small  element  of  surface  must  include: 

9*,  -  the  heat  flux  from  the  boundary  layer 

9m  -  heat  conduction  in  X  direction  within 

the  skin 

-  heat  conduction  in  Y  direction  within 
the  skin 

-  heat  loss  due  to  radiation  to  the  wind 
tunnel  (n)  and  to  the  model  interior  (m) 

The  term  ,  which  we  wish  to  determine,  should  be  the  dominant  term 
in  a  good  experiment. 


The  general  heat  balance  equation  is  written: 

y  c/yc/A  A  9*+*m)  -  9y,*r)t*« 

where  dW  is  the  element  of  mass  lying  beneath  the  surface  dAw  . 
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For  a  uniform  wall  thickness,  t,  and  density,  p m 


dW  =  dy  d  x  t 


Combining  all  elements  of  the  heat  balance  and  using  Taylor's 
expansion  and  Fourier's  law,  the  heat  balance  equation  becomes 


rm  c  -rr-  -  (  s~rr  *  *rrrx  / 


c/r  1 

^  (nC-rS) 


(i) 


2.  Estimate  of  Significance  of  Radiation 

It  is  possible  to  obtain  reasonable,  but  not  precise,  estimates 
of  the  radiation  heat  fluxes.  Radiation  losses  to  the  test  section  walls 
and  nozzle  blocks  must  be  considered;  however,  radiation  exchange  into 
the  interior  of  the  model  will  not  be  taken  in  account.  In  consideration 
of  the  extrapolation  technique  later  explained  and  the  uniform  temperature 
distribution  of  time  f"  =  0,  the  error  made  under  this  last  assumption  is 
negligible. 

In  order  to  estimate  the  heat  exchange  by  radiation,  we  assume 

the  model  walls  are  perpendicular  to  the  wind  tunnel  walls.  The  heat 

exchange  by  radiation  in  this  case  will  be  larger  than  that  corresponding 

to  the  actual  model.  For  this  case  the  maximum  radiant  heat  flux  for  the 

2 

hypothetical  surface  was  found  to  be  291  Btu/ft  -hr.  Since  the  smallest 
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initial  temperature  potential  was  170°F,  the  radiation  heat  transfer 

2  0 

coefficient  291/170  =  1.  71  Btu/ft  hr  F  was  at  most  4%,  (see  Figures  21 
through  26)  which  lies  within  the  error  range  expected.  Based  on  this 
calculation,  radiation  is  eliminated  from  further  consideration. 

3.  Estimate  of  Significance  of  Conduction 

As  was  stated  before,  a  uniform  temperature  distribution  along 
the  walls  of  the  bars  was  necessary  in  order  to  eliminate  heat-conduction 
within  the  skin  at  the  time  f  -  0.  Because  of  the  influence  of  the  free 
convection  and  of  the  non-symmetry  of  the  model  with  respect  to  each 
point  where  data  has  been  taken,  it  was  not  easy  to  obtain  the  desired 
uniform  temperature  distribution  during  the  heating  process.  The  largest 
temperature  difference  measured  at  the  start  of  the  run  within  a  cross- 
section  was  10°F.  On  this  basis,  the  largest  of  the  terms  &  Tw I  £  x 
and  d  Tw/  £  y  in  Equation  1  was  found  to  be  less  than  2.  0°F/ft  .  With 

2 

this  value  the  conduction  term  in  Equation  1  was  less  than  0.  008  Btu/ft  hr. 
Since  is  of  the  order  of  magnitude  10*  Btu/ft^hr,  the  conduction 

term  can  be  neglected. 

4.  Experimental  Technique 

With  the  help  of  Equation  1  with  the  conduction  and  radiation  terms 
dropped,  the  experimental  technique  can  be  discussed.  The  model  was 
heated  to  around  250°F  and  held  until  the  temperature  distribution  became 
uniform.  The  runs  were  made  during  a  typical  Minnesota  summer  day  during 
which  the  tunnel  temperature  averaged  80°F.  Hence  the  initial  temperature 
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potential  was  about  170°F,  When  the  model  temperatures  were  stable 

and  uniform,  the  model  was  inserted  into  the  wind  tunnel.  Model 

temperatures  decreased  with  time  as  a  result  of  the  convective  cooling 

process.  Each  term  in  Equation  1  is  time  dependent.  Initially  only  the 

storage  term  and  the  model  nozzle  radiation  terms  contribute.  As  time 

proceeds,  the  other  terms  increase  in  importance.  The  equilibrium 

situation  occurs  when  the  first  term  on  the  right  reduces  to  zero,  and  the 

model  is  at  the  adiabatic  wall  (T  )  condition. 

aw 

The  accuracy  of  the  method  depends  to  a  large  extent  on  the 
frequency  with  which  reliable  thermocouple  readings  are  recorded  near 
th <i  beginning  of  the  run.  Hence  measurements  were  taken  at  only  one 
cross-section  at  a  time,  with  temperature  readings  at  each  point  every 
1.4  seconds  for  the  cross-sections  B-B,  C-C,  B'-B',  C'-C'  and  every 
1.  2  seconds  for  the  cross-sections  A-A  and  A'-A'. 

Taking  data  early  in  a  run  implies  that  no  surface  heat  conduction 
corrections  are  necessary;  a  fact  which  greatly  simplifies  data  reduction. 
The  aerodynamic  heat  flux  is  then  obtained  for  the  isothermal  surface 
temperature  case,  which  is  the  most  basic  and  easily  specified  case  and 
the  only  case  for  which  most  heat  transfer  analyses  apply. 

5.  Data  Reduction 

Equation  1  is  now  written  in  the  form 

%/z  '  **  *  ~  ~ 


(2) 
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For  h  constant  and  T  fixed  by  the  prescribed  pressure  ratio, 

aw 

integrating  Equation  2  yields 


r  e 


/>  C' 


(3) 


where  C*  corresponds  to  the  elapsed  time  from  starting  when  wall 
temperature  was  Tw(0).  T&w  in  Equation  3  is  computed  using  the  ratio 
of  adiabatic  wall  temperature  to  stagnation  temperature,  defined  by 


r*  z 


r0 


This  ratio  is  measured  by  running  for  a  sufficient  length  of  time  (see 
Section  IIIC).  Thus  the  adiabatic  wall  temperature  T  can  be  calculated 


T 


aw 


r 


where  Tq  is  the  stagnation  temperature  measured  for  each  run.  Solving 
Equation  3  for  h 


h  - 


) 


(4) 


where  Tw  corresponds  to  the  temperature  at  time  T". 

Running  time  was  extended  for  approximately  20-30  seconds  after 
insertion  of  the  grid  assembly.  This  technique  allowed  a  time-temperature 
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history  to  be  recorded.  Then,  using  this  data,  the  heat  transfer 
coefficient  could  be  computed  at  any  time  from  Equation  4.  However, 
due  to  the  inherent  error  in  determining  the  elapsed  time  from  start, 
two  consecutive  data  points  were  used.  Applying  Equation  4  to  two 
different  times,  and  ^  and  reducing  we  get 


In  Equation  5  the  time  difference  appears.  This  was  very 

accurately  measured  by  the  Dymec  system. 

Seven  data  points  were  used  to  calculate  heat  transfer  coefficient 
values,  two  points  for  each  h,  yielding  six  different  coefficients.  These 
six  heat  transfer  coefficients  should  be  the  same  if  all  other  modes  of 
heat  transfer  are  absent.  However  due  to  conduction,  both  radially  and 
axia  uy.  the  coefficients  changed  with  time. 

Thus  it  was  decided  to  plot  these  calculated  coefficients  against 


a  dimensionless  temperature  ratio  Redefined  as 


'ah' 


where  T  is  simply  a  ^  metric  mean  of  the  two  successive 

'"Aver 

temperature  responses. 

Since  at  the  instant  of  insertion,  the  heat  transfer  mode  is 
entirely  convection  (temperature  distribution  uniform  about  periphery). 


I 

M 
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the  most  accurate  values  of  h  are  obtained  for  small  times,  i.  e.  ,  when 
Ry.  is  near  unity. 

By  plotting  the  heat  transfer  coefficient  h  versus  R  we  obtain  a 
history  of  the  convection  process.  (Sec  Figure  33).  The  data  points 
were  fitted  with  a  least  squa  re  polynomial  with  the  boundary  condition  of 
zero  slope  for  equal  to  unity.  This  extrapolation  technique  then 
afforded  a  heat  transfer  coefficient  independent  of  any  conduction  effects. 

This  method  of  analysis  proved  adequate  for  the  majority  of  runs. 
However,  for  some  runs  where  the  initial  point  was  inconsistent,  the 
time  temperature  histories  were  graphically  plotted  and  Equation  2  was 
applied  directly.  (See  Figure  34). 

The  extrapolated  values  or  hand  plotted  ones  are  then  the  heat 
transfer  coefficients  used  throughout  this  report. 


\ 
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FIGURE  33  TYPICAL  EXTRAPOLATION  TECHNIQUE  FOR  HEAT  TRANSFER  COEFFICIENT 
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